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Abstract 

With the eigenfunctional theory, we study a general interacting electron system, and give a rig- 
orous expression of its ground state energy which is composed of two parts, one part is contributed 
by the non-interacting electrons, and another one is represented by the correlation functions that 
are controlled by the electron correlation. Moreover, according to the rigorous expression of the 
ground state energy, an effective method beyond the local density approximation of the density 
functional theory is proposed. 
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The electron correlation plays an important role in description of many-electron systems. 
According to the electron correlation strength, the systems are approximately divided into 
two categories, one is called weakly correlated electron systems, and another one is called 
strongly correlated electron system [if]. The former can be approximately described by the 
independent-particle (quasiparticle) schemes [2|, yj, while the latter has to be represented by 
new schemes beyond the independent-particle description. 

The Hohenberg-Kohn-Sham (HKS) density functional theory^, 0] opens a new era in de- 
scription of the many-electron systems, and it and its generalizations [6|] are extensively ap- 
plied in physics and chemistry. However, the rigorous expression of the exchange-correlation 
energy in the HKS density functional theory is unknown, and the local density approxima- 
tion (LDA) of the Kohn-Sham scheme plays a central role in calculating the ground state 
energy of the systems. In general, for the weakly correlated systems, the result of the LDA 
is accurately agreement with experiment data, but it is unreliable for the strongly correlated 
systems because of the strong electron correlation. 

In contrast to the Kohn-Sham scheme of the density functional theory, we study a gen- 
eral quantum many-electron system with the eigenfunctional theory 3] , and give a rigorous 
expression of the ground state energy, especially the correlation energy. Moreover, according 
to this rigorous form of the ground state energy, the correlation energy can be well defined, 
and a more effective method beyond the LDA is proposed. 

The Hamiltonian of a many-electron system can be generally written down, 

~ p 2 r 

H = H + — d 3 xd 3 x'V(x - x')p(x)p(x') (1) 

where H = J dxip^x) (j^ + U(x)^j ip a (x) is the Hamiltonian of the non-interacting 
electrons, ip a (x) (ip a (x)) are the creation (annihilation) operators of the electrons with spin 
label cr at the coordinate x, = '^2 cr ip^ a (x)ijj a (x) is the density operator of the elec- 

trons, and U(x) is an external potential. The last term represents the electron Coulomb 
interaction, and it induces the electron correlation effect which makes the problem hard be 
treated. In usual mean field theory 0,3, the (effective) Coulomb interaction is treated as 
a perturbation parameter, while in the density functional theory, it is incorporated into the 
exchange-correlation energy. However, if the electron correlation is strong, in which case the 
results obtained by the mean field theory and the LDA of the density functional theory are 
unreliable, and the electron Coulomb interaction has to be accurately treated. 



With the eigenfunctional theory, the partition function of the system reads Q], 

Difj*DipDpD(j)e* s 
S = [ dxdtip* a (x,t)M(x,t)^ a (x,t) +W\p,4>] (2) 

where W[p, 0] = / dxdt(p(x, t)p(x, t) — y / dtd 3 xd 3 x'V(x — x')p(x,t)p(x',t), M(x,t) = iHd t + 

p — 2 U (x) — (p(x, t) is the electron propagator operator in the new Hilbert space, and p 

is the chemical potential. The Lagrange multiplier field <fi(x, t) is introduced to decouple the 
electron interaction, in the meanwhile the original Hilbert space is mapped into a new Hilbert 
space, in which the electrons are non-interacting, and moving in a fluctuating potential 
produced by 6{x, t). 

u n 

In general, the eigenequation of the electron propagator operator reads [71. 1 101]. 

M(x,t)* aku (x,t; [4>]) = E aku [<j>]* aku (x,t; [4>]) (3) 

and the eigenvalue can be obtained by the Hellmann-Feynman theorem, 

E akuj [4>] = Tuo - E k - E ak [(j)} 
Z ak [4>] = dX dtd 3 x<t)(x,t)\y akuj {x,t;[\<t)))\ 2 



o 



where uo is the frequency, H ak [4>] is the self-energy of the electrons in the new Hilbert space, 
and E k is the eigenvalue of the non-interacting Hamiltonian of the electrons, 

H ip ak (x) = (E k + p)iJ ak (x) 

where H = + U(x), and k labels a set of quantum number representing the states of 
the non-interacting electrons. According to the expression of the eigenvalue E akw [<p], the 
eigenfunctional can be generally written down, 

<W*, *5 M) = -^^(^e^^M^^Wi (4) 
VT 



where T — > oo is the time scale of the system, and the phase field Q ak (x,t] [0]) satisfies the 
eikonal-like equation with the condition Q ak (x,t; [<fi = 0]) = 0, 

(m - £ - ^IpH^m-p) Qak(x,t; [0]) - ggf^ ^ ])]2 = 0( M ) (5) 
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For the homogeneous case, this equation can be easily solved after neglecting the non 



linear term which in general is a small quant ity[ll|]. It is worthily noted that the difference 
between the eigenfunctional ^ a ku>( x i t] [0]) of the electrons in the new Hilbert space and 
the wave function ifj ak (x) of the non-interacting electrons is the functional eQ"*>( x > t '>l<l > ]) but a 
pure phase factor e -l( - a, ~ So " fc [™*, thus the physical meaning of the phase field Q a f.(x,t; [0]) is 
very clear that the electron correlation effect can be completely represented by the phase 
field, i.e., the phase field is a correlation parameter of the electrons. The quantum phase 
transition of the system from the Landau Fermi liquid to the non-Fermi (Luttinger) liquid 
can be described by the phase field. 

In terms of the eigenfunctionals ^ a ktj( x , t; [0]), the second quantization representation of 
the electrons in the new Hilbert space can be written down, 

VvO,t) = y^^ aku} (x,t; [</)]) c akul 

ku 

fl(x,t) = X>; to 0M;M)4fa, (6) 

ku 

where c a fc w (stku) * s the annihilation (creation) operator of the electrons with the spin index 
and the quantum number k and oj. Using the orthogonality and the completeness of the 
eigenfunctionals ^/ a ku(x, t; [0]), it can be easily proved that the above electron operators 
satisfy the standard electron anti-commutation relations. These expressions play central 
role in calculating the ground state energy and a variety of correlation functions of the 
system. 

After integrating out the electron fields, and using the mathematical formula, 

f 1 1 

Tr In (A + b) = Trln ( A) + Tr / dXB— - 

A + XB 



o 



the partition function can be written down 
Z = J DpD(j)e^ s[p ^ ] 

S[p,(j)} = -iTrln(M ^ +W[p,<j)] + iJ2 J Q dX j 



dtd 3 X(f)(x, t)G a (x, t; x, t; [X<p]) (7) 



where Mo = ihdt + p — — U (x) is the non-interacting electron propagator operator. The 
electron Green function G a (x, t; x' , t'; [0]) in the new Hilbert space can be written down, 

C ( r f- r > f >. U]) - V gg*!fe!i t^)^^' 
u a {x,i,x ,i , m) - ^ h _ „ _ T rn 



and the electron Green function in the original Hilbert space is that, G a (x,t; x' ,t') =< 
G a (x,t;x',t'; [(f)}) > p<j) , where < A[<j>] > p4> = \ J DpD<j>e^ s ^A[<j>]. The action S\p,<j>] is used 
to calculate the functional average over the auxiliary fields <f>(x,t) and p(x,t). 

According to the expressions of the electron operators, the ground state energy E g =< 
H(t) > of the system reads, 

E g = T + Vhf + AT + E c (8) 

where T = 2 J2k n k (^k + P) is the kinetic energy of the non-interacting electrons, and 
n k — @(—Ek) is the occupation number of the k-th state of the non-interacting electrons. 



The term Vhf is the usual Hartree-Fock potential energy 
e 2 f 

Vrf = y / d 3 xd 3 x'V(x - x') p (x) p (x') 



12 



e 2 



akk' 



— ^2n k n k , I d 3 xd 3 x'V(x - x')^l k (x)ifj ak ,(x)ifj* ak ,(x')ifj ak (x') (9) 



where p (x) = J2ak a ki. x )\ 2 is the density of the non-interacting electrons. The first 
two terms in the expression of the ground state energy is similar to that in the LDA of the 
density functional theory, which can be rigorously represented by the single-particle wave 
functions. 

The term AT is a modification of the kinetic energy of the electrons by the Coulomb 
interaction |l3|. 

AT = J2j d 3 xp ah (x)g ak (x, [p ]) (10) 

ak 

where p ak (x) = n k \ip ak (x)\ 2 , and gak(x,[p ]) =< \e Qak{x ^ ]) \ 2 [ihd t Qak(x,t][(j)\) - 
<j>(x,t)] > fX j > is the functional of the density p (x), and independent of the time coordinate t 
because of the time translation symmetry of the system. For a homogeneous electron system, 
the function g ak (x, [p ]) = g ak [p Q ] is equivalent to the self-energy under the random-phase 
approximation, and it is a functional of the density p {x). Thus, the term AT can be treated 
approximately under the LDA. 



The correlation energy E c = E Hc + E Fc can be written down 



e f 

E Hc = y / d 3 xd 3 x'V(x - x')p(x,x') 
e 2 r 

E Fc = — d 3 xd 3 x'V(x-x')r(x,x / ) (11) 
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where the functions p(x,x') and F(x,x') are defined as, 
p{x,x') = ^2 n k n k/ \ifj ak {x)\ 2 \if) pk (x')\ 2 

a/3kk' 

r(x,x') = ^n k n kl ^l k (x)^ akl {x)^l kl {x')^ ak {x) 

akk' 

Obviously, Eh c is the Hartree-like correlation energy, while Ep c represents the Fock-like 
correlation energy. The former shows the local behavior, and the latter has the non-local 
property. It would be demonstrated that the main errors of the LDA of the density functional 
theory originates from the correlation energy Ep c term. 

For a homogeneous interacting electron system, according to the experience in studying 
the one-dimensional interacting electron systems [17]], the electron correlation strength is 
controlled by the imaginary part of the phase field Q a k{x,t; [0]), and it can be shown that 
< | e Qrt(a,t;M)|2| e g pfc («',tiM)|2 >p(j)= e F aPkk ,{x-x';iP }) } w h e re the function F a p kk ,(x - x'\ [p ]) is a 

small quantity, thus it is reasons to ta,e the app—ionQ ~ 

1 + F. Therefore, 

the function p(x, x') can be approximately represented by the density function, 

p(x, X 1 ) ~ ^ F ^kk'( X ~ X 'i [Po\)Pak( X )p/3k'( X ') ( 12 ) 
a(3kk' 

where F a ^k'( x ~ x '] [p ]) is the functional of the non-interacting electron density p Q (x). Thus, 
the correlation energy Eh c can be treated approximately under the LDA. 

It can be shown that < e^k^iMJ+^^AMJ+^^AMJ+^v^iM) > p0 = ^^-^M for 
the homogeneous case, while the function P akk i(x — x'; [p ]) is not a small quantity for the 
systems with the strong electron correlation. For example, for one-dimensional interacting 



electron gas, the function e p ° kk '( x x '^ol) has a power-law asymptotic behavior 19], 

e Pl°,(*-*.';[p o] ) 1 (13) 

\x-x'\~< v ' 

where 7 > 1 is a dimensionless coupling constant which depends upon the electron interac- 
tion strength and the Fermi velocity. In this case, the function T(x,x') shows completely 
different behavior from that of the function p(x,x'), and the correlation energy Ep c cannot 
be treated approximately under the LDA. In general, it is expected that the correlation en- 
ergy Ep c is important for the systems with the strong electron correlation, and it dominates 
the correlation energy of the systems, thus a method beyond the LDA is needed. 



The above expression of the ground state energy is obviously different from that of the 
density functional theory in the LDA (Kohn-Sham scheme) where the ground state energy 
is represented by the quantity of the non-interacting electrons. The ground state energy 
in equation (jSj) is composed of two parts, one part is corresponding to the non-interacting 
electrons, and another one is represented by the correlation functions that are controlled by 
the electron correlation which describes the many-electron effect. In principle, the expression 
of the ground state energy is independent of the choice of the Hamiltonian Ho of the non- 
interacting electrons. The different choice of the Hamiltonian H can give different non- 
interacting part, while for the correlation part, it will only change the action S[p, 4>] and the 
wave function if) ak (x) of the non- interacting electrons. 

Based upon the above discussions, a more effective method for calculating the ground 
state energy E g is that, by suitably choosing the non-interacting Hamiltonian Hq, the terms 
Tq, VhF) AT and Eh c can be calculated with usual LDA scheme, then the term Ep c can be 
calculated by the quantum Monte Carlo method (QMCM) 2(| where the action S[p, <fi] can 
be obtained by solving the equation (J5j) of the phase field. By more accurately treating the 
correlation energy Ep c , the results of usual LDA will be heavily modified for the systems 
with strong electron correlation. 

On the other hand, for a suitable LDA, if the contribution of the correlation energy Ep c 
to the ground state energy is small, this LDA is reliable. If its contribution is important, 
a new LDA has to be chosen to make its contribution be as small as possible. Of course, 
corresponding to the local spin density approximation (LSDA) 16j of the density functional 
theory, it needs to introduce the Lagrange multiplier fields (p a (x,t), in which way there 
are four auxiliary fields (p a (x,t) and p a (x,t), while the ground state energy has the similar 
expression as that in equation (jSJ. 

In summary, with the eigenfunctional theory, we have studied a general interacting elec- 
tron system, and given a rigorous expression of its ground state energy which is composed of 
two parts, one part is contributed by the non-interacting electrons, and another one is repre- 
sented by the correlation functions that are controlled by the electron correlation. Moreover, 
according to the rigorous expression of the ground state energy, an effective method beyond 
LDA may be the LDA plus the QMCM. The LDA is used to treat the non-interacting 
electron part, while the correlation part is treated by the QMCM. 

This work is supported by the NSF of China, No.90403015 and No. 10774188. 
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